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ABSTRACT 


Current meter data were collected from moorings on the 800 and 1800 m isobaths 
on the continental slope off Point Sur, California at 100, 350, and 500 m depth and at 
1000 m depth on the 1800 m isobath. Complex demodulation and spectral techniques 
were used to study the distribution of tidal energy at the diurnal and semidiurnal fre- 
quencies, and to study the size, rotation and orientation of the tidal current ellipses. At 
moonng P2 (800 m) the observed diurnal energy increased with depth for each of the 17 
months of data analyzed, and the semidiurnal energy likewise increased with depth for 
14 of the 17 months. The observations were studied using the theory of bottom trapped 
waves in a rotating stratified fluid (Rhines, 1970). Bottom trapping is expected at the 
diurnal and semidiurnal frequencies whenever the critical trapping frequency 
w, = Nsinasing exceeds 1/24 and 1/12 hours"! respectively, where N = the Brunt- Váisálá 
frequency, a the bottom slope, and $ is a measure of the angle the wavenumber vector 
makes with the gradient of topography. Theoretical energy decay curves matched the 
observed energy levels at P2 very well, indicating that the enhancement of energy with 
depth was likely due to the presence of bottom trapped waves. No bottom trapping was 
observed at P3 (1800 m) because the deepest current meters were positioned too far 
above the sea floor to detect bottom trapped energy, and because the density gradients 


near the sea floor were weaker than at site P2 and kept c, below the tidal frequencies. 
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l. INTRODUCTION 


A. BACKGROUND 

The Naval Postgraduate School maintained an array of subsurface current meter 
moorings off Point Sur, California from 1959 to 1990. The moorings were positioned 
on the continental slope and collected data used in the study of the California Current 
System (Figure 1). The initial spectral analysis of data from mooring P2 indicated that 
the kinetic energv at diurnal tidal frequencies increased with depth (Figure 2). This was 
an unusual phenomena because vertical energy structures at tidal frequencies are usually 
associated with the baroclinic modes of free internal waves. At the latitude of Point Sur 
(36° 20’ N) the inertial frequency is approximately 1.1 cpd and free internal waves can- 
not exist at diurnal (subinertial) frequencies. Baroclinic waves can exist at subinertial 
frequencies only if the waves are trapped along a wall or slope. 

The tides off most of the California coast are a mixture of diurnal and semidiurnal 
components. The principal lunar constituent, M, is the strongest tidal element followed 
by the A, O, and S, constituents. Barotropic diurnal and semidiurnal tides both propa- 
gate northward along the coast. The semidiurnal tide is a combination of Kelvin, 
Poincare and forced waves. The barotropic semidiurnal tidal current ellipses, represent- 
ing horizontal current amplitudes and directions, are generally oriented parallel to the 
coastline (Munk et al., 1970). The barotropic diurnal tide is a combined Kelvin and shelf 
wave with current ellipses which can be either parallel or perpendicular to topography 
(Noble el al., 1987). Several studies of California tides have found the primary direction 
of rotation for diurnal and semidiurnal current ellipses to be counterclockwise. 

Baroclinic tides are generated when semidiurnal surface tides force isopycnal sur- 


faces over topographic features, such as the continental shelf break (Wunsh, 1975). 
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Figure 1. Point Sur Moorings P2 and P3: The location of Point Sur on the 
California coast is shown in the inset. 


While the barotropic semidiurnal tide is usuallv aligned with topography, internal tides 
often become the dominant cross-slope current component along the California coast 
(Noble et al., 1987). 

Variations in the vertical distribution of tidal energy have been observed at 
semidiurnal frequencies. During the Coastal Ocean Dvnamics Experiment (CODE), the 
barochnic semidiurnal internal tide. was larger than the barotropic tide. Changes in 
semidiurnal kinetic energy at different depths were linked to surface winds and changing 
stratification (Rosenfeld, 1990). The semidiurnal internal tides observed in CODE had 
some bottom energy intensification which was attributed to the horizontal velocity 


structure of the first dynamic mode. 
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Figure 2. Spectral Analysis of Energy at Mooring P2 During July 1989: The 
curves represent energy at 100 m, 350 m and 500m. 


There are few published studies which discuss similar bottom intensification of 


subinertial currents. Bottom trapped topographic Rossby waves were observed by 


Thompson and Luyten (1976) and were discussed using bottom trapped wave theory 
developed by Rhines (1970). These were low frequency waves, however similar theory 
was applied to bottom trapped waves at tidal frequencies and used to study the tidal 


currents and vertical energy distribution off Point Sur. 


B. BOTTOM TRAPPED WAVES IN A ROTATING STRATIFIED FLUID 

This discussion follows the development presented by Rhines (1970) for bottom 
trapped Waves in a uniformly stratified, rotating fluid bounded by a single plane wall 
inclined everywhere at angle « to the horizontal. A right handed orthonormal coordi- 
nate system 1s oriented on the slope such that the x, and x, axis are rotated at an angle 
$ from the alongslope and upslope directions respectively and x, 1s perpendicular to the 
bottom slope (Figure 3). The corresponding velocities along x, x, and x; are designated 
as u, , u, and t. 

The governing equations are the linearized, non-diffusive Boussinesq equations used 


in the form: 





à Vp „A 

20 cen (1) 
Cp’ N? A 
z =g Y (2) 


where p = pressure, Q = angular rotation of the earth, g = gravitational acceleration, 


j = a vertical unit vector, 22 = ff = f(sinasing, Sinacosý, cosa) , f= 
- —g 0 
Coriolis parameter, the Brunt-Váisálà frequency squared 2 Nº = ENS , and a per- 


turbation density is used such that p = p, + p(z) + p’(x, y, z, 1). 
Inviscid boundary conditions are chosen so that uj=0 at x; =0. If u =0 every- 


where in the water column the motion is nondivergent in planes perpendicular to the x; 


Ó O : : 
axis and the continuity equation = F E — 0 applies. The wave propagates in the 
l 2 





Figure 3. Coordinates for Bottom Trapped Waves on a Sloping Bottom: The ele- 
vation from the horizontal is given by a. 
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x, direction with = O. To balance the continuity eguation ET must also equal 
l 2 


zero. When u,#0 and X-- O, free internal waves can exist. If u, — 0 and É 
2 2 


waves are trapped against the wall (bottom) in a plane parallel to the boundary surface 
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with particle motion limited to the x, direction. 
Substituting Uu, u, = u, — 0, and the appropriate form of 2Q into Equations (1) and 


(2) leads to a set of three equations where for the x, direction, Equation (1) becomes 
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In the x, direction, Equation (1) becomes 
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For the x, direction, the result is 
1 Cp , 
— fu,( sin a cos $) 2 — poe m go'( cosa) (5) 


The rate of change of the density perturbation, equation (2), similarly reduces to 
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The combination of equations (3) and (6) results in an expression for the frequency a 


where: 
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a=+Nsinasiné (7) 


Substituting (7) into (6), p' can be determined as 


Solutions to the wave eguations are found using the form 


p= poe * + lk3)x, + (ky Xx, — ot) (9) 
where x and Á, are real. x is defined as a trapping coefficient with units of wavenumber. 
The alongslope or horizontal ( k, ) and vertical ( k, ) wavenumbers in this solution lead 
to spatial dependence in the x, and x, directions. Using equations (8) and (9) in (4) and 
(5), the spatial dependence of the wave 1s: 
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N , est 
where S = a When the trapping coefficient x is positive the wave motion is expo- 
nential and trapped near the bottom. From equation (9), wave motion decreases due to 


the effect of trapping by 
CR (12) 


The magnitude of trapping is dependent upon o and q. When particle motion is di- 
rectly up the slope (¢ = 90°) trapping 1s maximized. When particle motion is along the 
isobaths (@ = 0°), no trapping occurs. 

Rhines discusses four different geometries where trapping occurs: (1) a ~ 1 for steep 


slopes, (2) a ~ S7! for more gradual slopes such as those on the continental slope, (3) 


a< $! for very gradual slopes such as those on the continental rise, and (4) 
T : i : ; 
> —a<Scos @ for near vertical slopes. Consider the case appropriate for the conti- 
nental slope off the California coast where a ~ S- and sino ~ O(1). Evaluating 


equation (7) using these relationships leads to the following: 


o = LAN Sing SI d ~ N sin(S7') sin 6 ~ ve 
í 
of (13) 
= can be expressed as: 
" sf. ane sind 
= = E RS CT o = 5 (19) 
K Cos Sms SILES” MO 
K; 
The ratio of To reduces to 
ke SEI — sint(S7?) sin’é J 
A = a o (15) 
ky cos(S )J+Ssin(S )cos o 
The relationships now established with a~S-' are o-f, ES, and 


k, 
k 
— < S. Parallel to the slope Coriolis is balanced by buoyancy as shown in equation 


A 

(15). Coriolis traps waves against the wall (bottom) but the frequency of these trapped 
Waves, determined by. equation (7), 1s not dependent on f. Bottom trapped internal 
waves are similar to internal Kelvin waves because they are not bound bv a lower fre- 
quency limit of. f and trapped waves can exist and propagate at subinertial frequencies. 
The high frequency limit of bottom trapped waves is o and is a function of a and N. 


Maximum trapping will always occur when ¢ = 90° and the particle motion is directly 


up the slope. 


The exponential decay of trapped waves away from the bottom is determined bv 
trapping coeíficients. To determine x, bottom trapped waves were assumed to have a 


dispersion relationship similar to baroclinic Kelvin waves (Rhines, 1970) where 


O. 
DEE Nk 
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with w the frequency of a wave, k a horizontal wave number, n the mode number and 
H the water depth. The use of this dispersion relationship represents an approximation 
because baroclinic Kelvin waves assume a flat ocean bottom. The horizontal wave 


number can then determined using 


m 


k= wh 2 2 
(A = p) 


When À = 90° waves are propagating along the slope in the x, direction and particle 
motion is up and down the slope in the x, direction. Trapping coefficients are defined 
as kK = k,S from equation (14) and can be calculated using 


E (16) 


f (N? — w?) 


Waves are trapped closer to the sea floor and decay more rapidly away from the 
bottom when x increases. Trapping coefficients increase when (1) the buoyancy fre- 
quency N increases, (2) depth H decreases, (3) higher frequencies w are considered (in 
the tidal range of 1.0 to 2.0 cpd), (4) f decreases (lower latitudes) and (5) for higher 
modes n. Trapping 1s weakest and waves decay more gradually away from the bottom 


in the first baroclinic mode. No trapping occurs for the barotropic case when n = 0. 


IL DATA COLLECTION 


A. CURRENT METER DATA 
Current data were collected at two sites off Point Sur, California (Table 1) located 


over the continental slope (x = .06). 


Table 1. POINT SUR MOORING LOCATIONS 


MO l EE 
36-20.0 N 122-10.2 W 800 m 
36-20.0 N 122-27.6 W 1800 m 







Mooring P2, on the $00 m isobath, was equipped with three Aanderaa vector averaging 
current meters (RCMS8) at depths of 100 m, 350 m and 500 m. Mooring P3 was situated 
52 kin offshore in 1800 m with four Aanderaa current meters at depths of 100 m, 350 
m, 500 m and 1000 m (Figure 4). 

Vector averaged current velocities, current directions, temperature, conductivitv and 
pressure measurements were recorded every 30 minutes in a solid state Data Storage 
Unit (DSU). The DSU had a storage capacity of 65530 10-bit words and could record 
seven months of data with the sampling interval used. Sampling times were regulated 
by an internal quartz clock with an accuracy of +2 seconds per dav. Several deploy- 
ments were necessitated by the seven month data storage limit (Figure 5). P2 was re- 
covered and redeployed m May, August and December 1989 and May 1990. The three 
current meters used on P2 provided nearly continuous data from Mav 1989 through 
October 1990. Mooring P3 was deployed in December 1989 and again in May 1990. 
Good data was obtained at four levels from December 1989 to April 1990. After April, 


the data from the current meters at 500 m and 1000 m was limited by instrument failures. 
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Figure 4. Current Meter Configuration at Moorings P2 and P3 


Velocities were measured using a shrouded paddlewheel at the top of the recording 
unit. The paddlewheel was magnetically coupled to an electronic counter and velocities, 
in units of cm sec”?, were determined by the number of revolutions the paddlewhcel 
made during each sampling interval. The sensor had a range of 2 - 250 cm sec”! with 
an accuracy of + 1.0 cm sec”! or + 2 percent of the current velocity. The paddlewheel 
Stalled at velocities below a threshold of 2.0 cm sec”! (Aanderaa, 1990). The cosine re- 


sponse was acceptable within + 27° from the vertical. Both the pressure sensors and the 
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Figure 5. Current Meter Deployments at Moorings P2 and P3: The shaded bars 
indicate the intervals where good data was obtained at each depth. 


mooring dynamics program in use at NPS indicated a maximum moonng tilt of 11.7“ 
which was well within acceptable limits. 

The current meters were oriented into the current by a large vane with an area of 
2425 cm?. The direction sensor was a magnetic compass with a potentiometer located 
at the bottom of the recording unit. The compass was oil damped and required 3 to 5 
seconds to resolve changes of 90° or more. The resolution of the compass was 0.35° 
with an accuracy of + 5° for velocities 5 - 100 cm sec-' and + 7.5 degrees for velocities 
of 2.5 - 5.0 cm sec"! and 100 - 200 cm sec-!. The compass could function properly with 
a maximum tilt angle of 12° from the vertical. Prior to each deployment the combined 
magnetic deviation and variation of each unit was determined on a surveyed test bench 
at the Naval Postgraduate School and the necessary corrections were applied via 


postprocessing software. 
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B. HYDROGRAPHIC DATA 

While the current meters Were deploved, several oceanographic research cruises were 
conducted in the region by the Naval Postgraduate School. These vovages consisted of 
California Undercurrent (CU) and Student Cruises (ST) which collected CTD data at 
intervals of approximately 10 km along the Point Sur Transect (POST). Data were col- 
lected near the moorings from the surface to within 50 - 150 m of the bottom using a 
Neil Brown Instruments System (NBIS) Mark IIIB CTD and are described in detail in 
Tisch (1990). Using the data from the cruises (Table 2), the various characteristics of 
the waters off Point Sur were examined. Although the station locations and spacing 
varied slightly with each cruise, the CTD data necessary for interpretation of the current 


meter results were comparable. 


Table 2. HYDROGRAPHIC CRUISES OCCURRING DURING 
CURRENT METER DEPLOYMENTS 
















Ill. METHODS OF ANALYSIS 


A. SMOOTHING AND FILTERING 

Upon recovery, the current meter DSUs were downloaded into a computer. The raw 
data files were converted standard scientific units using the sensor calibration equations 
and stored as formatted ASCII files. All data records were plotted to facilitate visual 
inspections. Velocity records were truncated where instrument failures or malfunctions 
were evident. Spurious values, such as an occasional outlier, were edited or deleted 
manually. 

All data were band pass filtered prior to further analysis to isolate the frequencies 
of interest. The filtering and interpolation methods used were originally developed by 
Oregon State University (Denbo et al., 1984). Postprocessed records, containing data 
at 30 minute intervals, were first low passed filtered (LP) using a Cosine-Lanczos filter. 
Such a filter minimizes unwanted ringing or oscillations caused by the Gibbs phenomena 
while still providing a sharp roll off at the desired cutoff frequency. A centered 25 point 
segment of the records was used to produce each LP filtered data point. The LP filter 
had a half power period of 2.9 hours (or a frequency of 8.4 cpd). The filter allowed any 
signal with a period exceeding 3.6 hours to pass through unattenuated. Signals with a 
period of less than 2.0 hours were removed entirely. After the LP filtering, records were 
interpolated to transform 30 minute data records to 60 minute intervals with each data 
point faling on an even hour. The interpolation routine utilized Lagrangian 
polynomials to interpolate any unevenly spaced data points (due to missing individual 
data points or outliers) to specific time intervals (Gerald and Wheatly, 1989). After in- 
terpolation, the data was again low passed filtered (LLP). The second Cosine-Lanczos 


filter had a half power period of 46.59 hours (or a frequency of 0.515 cpd) and used a 
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centered 121 point segment to produce one filtered data point. This filter was designed 
to remove diurnal and shorter period tidal energy and any inertial frequency energy 
present in the signal. The LLP records were subtracted from the LP records to produce 
the band passed form used in the tidal analysis. A sample from 100 m at P2 (Figure 6) 


shows the results of the filtering procedures. 


B. TIDAL COMPONENTS 
1. Harmonic Analysis 

To understand the nature of the tidal currents, an analysis of the various tidal 
components was desired. Traditional harmonic tidal analvsis, or the technique of linear 
combinations, requires a complete set of observations over a standard length period. 
Tidal heights or motions can be expressed as a sum of harmonic terms, each determined 
by a known, fixed frequency. A single estimate of the amplitude and phase are then 
extracted through the analvsis of observations. Observational periods are normally 
multiples of synodic periods involving the most important tidal constituents. The pre- 
ferred record length for a complete tidal analvsis is 369 days. This allows extraction of 
the amplitude and phase of at least 64 of the largest tidal constituents (Dronkers, 1964). 
An adequate description of the major constituents can be accomplished using a much 
shorter observation of 29 days. A 14 day time series is adequate to provide the resol- 
ution of the principal diurnal constituents, K, and O,, and a 15 day time series will re- 
solve the principal semidiurnal constituents M, and S, Some of the longest tidal 
constituents require a full 19 vear record for accurate resolution (Schureman, 1958). 

The current records collected off Point Sur vary in length and are not ideally 
suited for the method of linear combinations. The traditional method could provide only 
a few independent estimates of the tidal constituents for each month of recorded current 
observations. Complex demodulation, a least squares analvsis technique, was used in- 


stead to accomplish the harmonic analysis. The primary advantage realized with this 
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method is that records of any length can be analvzed. Complex demodulation utilizes 
a set of normal equations which are solved for unknown amplitudes and phase corre- 
sponding to n frequency constituents. The number of unknowns solved for in the set 
of equations is 21 + I. The number of equations used in each set is determined by the 
number of observations, ri, utilized for each estimate and the solutions are determined 


using matrix algebra. Schematically, this process an be represented as 


[observed velocities] = [known tidal constants][unknown empirical constants] 


The least squares method minimizes the difference between an observed current velocity 
and the sum of the solutions corresponding to that observation. The difference, or error, 
can be attributed to unresolved tidal constituents, inertial currents or noise in the signal. 
The extracted solutions reflect only those current variations with a coherent phase and 
no assumptions are made about data outside the interval upon which the solutions are 
fit (Pugh, 1986). 

The analysis of diurnal and semidiurnal constituents was conducted starting 
every six hours over 48 hour segments and produced a large number of independent es- 
timates for each current record. The 48 hour period did not allow adequate separation 
of the O, and K diurnal and the M, and S, semidiurnal constituents which had prevented 
by synodic periods of 13.661 and 14.765 days respectivelv. Therefore the major udal 
cogsutuents (O, K, P, li Mi dy MS No and K,) were generalized as diurnal ( 
w “1 cycle per day) or semidiurnal (w ~ 2 cpd). Longer period constituents, such as the 
Lunar Long Period Mm and MSf, were removed by low pass filtering. Harmonics such 
as the Lunar Terdiurnal M, or Solar Terdiurnal S, were not considered due to the rela- 
tivelv small amplitude of these constituents when compared to the primary diurnal and 


semidiurnal constituents. Higher order harmonics Were removed by the initial band pass 


filtering. Inclusion of the inertial frequency (o = 1.1 cpd) as a known frequency pro- 
vided a better fit of the least squares solutions and reduced error. 
The complex demodulation program calculated the coefficients using a least 


squares approximation to the set of xı equations with 
uorv=A,+ A, sin w(t — t) -- 44 cos w(t — 1,)... (17) 


where o is the angular frequency of the tidal constituent, z is the time of each hourly 
record and /, is the time of the first data point of an entire current meter record. The 
number of terms can be expanded to A,,., to include all desired constituents n. For u 
velocities, data was placed in a single row matrix U. The coefficients of equation (17) 
are placed in matrix X. The following relationship was used to obtain the solutions 


to the set of equations represented by equation (17): 
(UR) A ) = A 


Similar calculations were done for the v velocities to create matrix AV. Once the coef- 
ficients A were determined (the values in matrices AU or AV), the amplitude and phase 


of each frequency component was calculated using 


Amp, = A) + d 


0 = aretan Lane 
An 


With these values, the current velocities can now be defined as 


u or v — Amp, sin(o(t — 1) 4- 0) 


Diurnal and semidiurnal amplitude and phase were extracted from each current 


meter record using the complex demodulation technique. An example is provided in 
mee 7. 


2. Construction of Current Ellipses. 


The amplitude and phase of clockwise (O, 0.) and anticlockwise (Q, 0,.) ro- 
tating vectors were computed using monthly averages of the component amplitude and 


phase vectors (U, 0,) and (V, 0,) obtained from complex demodulation using standard 
equations (Pugh, 1986): 


Q. 


o 2 + 
— [U + W — 2UV sin(0,— 6) [2 

] 72 2 r i + 
Ou = LU" + VÍ + 2UV sin(6, — 6,) |? 


9 U sin 0, -- V cos 0, 
SE 0, — V sin 6, 


a U sin 0, 4- V cos 0, 
= arctan A 
ge —U cos 0, — V sin 0, 
Using the rotary components, the dimensions of a current ellipse may be defined as 


Semi-major Axis = Q, + Qac 


Semi-minor Axis = |Q, — 0. 


Orientation of Semi-major Axis (clockwise from North) = bi = > (AO) 


Current rotation is clockwise if O, > Q,, and anticlockwise if Q,, — O.. Sample ellipses 


are shown for mooring P2 during December 1989 (Figure 8). 
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Figure 7. Complex Demodulation of the Diurnal Frequency at 500 m on Mooring 
P2: The first two time series plotted are amplitudes in cm sec”! and the 
next two curves are the phases of the u and v diurnal velocity compo- 
nents in degrees relative to the beginning of each record. 
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TIDAL CURRENT ELLIPSE 
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Figure 8. Tidal Current Ellipses: Current ellipses are shown for diurnal and 


semidiurnal tidal currents at P2 for December 1989. North is towards 
the top of the page. 
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Current ellipses from different moorings and varying depths were evaluated by 
visual comparison of the ellipse orientation, size, eccentricity and rotation and by 
quantitative comparison of the energy represented by each ellipse. Energv was com- 
puted in units of cm? sec? by the addition of the squared dimensions of the semi-major 


and semi-minor axis. 


C. SPECTRAL ENERGY ANALYSIS 

To evaluate the energy present in the tidal currents, a program which computes the 
spectral characteristics of a time series was employed. The program, VCSPC3, was used 
to calculate the kinetic energy (KE) spectra as the total KE spectra for both the u and 
v velocities of a current meter record over one month periods. The program utilized a 
Cooley-Tukev fast fourier transform (FFT). One month records were divided into five 
sequential, equal length segments. A separate spectrum was computed for each segment 
and the five segments were then averaged to produce the total energy spectrum. A 
parabolic data window was used in the program to eliminate discontinuities at the be- 
ginning and end of records created because record lengths were not exact multiples of 
the periods of each wave present in the signal (Bendat and Piersol, 1986). The effect of 
the data window was to improve the resolution of the narrow spectral peaks in the tidal 
frequencies by reducing spectral side lobe leakage and associated distortion. No addi- 
tional smoothing of the spectral estimates was done. The program calculated 95% 
confidence intervals using a Chi squared distribution. 

Program VCSPC3 provided both graphical output and data listings. The plots show 
kinetic energy density in units of cm? sec? cpd^! versus frequency in cpd. Confidence 
intervals are indicated bv the separation of parallel horizontal lines above the horizontal 
axis. Distinct peaks were clearly seen at both the diurnal and semidiurnal frequencies 


of these plots (Figure 9). 
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Figure 9. Spectral Energy at P2: A representative plot from the spectral analysis 
program is shown for a one month interval. The 95% confidence inter- 
vals is plotted below the curve. 
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The listings show spectral density at the diurnal and semidiurnal frequencies with a 


frequency resolution of 0.05 cpd and the associated 95 percent confidence intervals. 


D. HYDROGRAPHIC DATA 

A series of temperature (T), salinity (S), density (o), geostrophic height (©), 
geostrophic velocity (W), and spiciness (7), values were plotted in contoured vertical 
sections Which passed through the locations of moorings P2 and P3 (figures not shown). 
A level of no motion for the geostrophic velocity calculations was chosen at 1000 m 
except where bottom depths were less than 1000 m when the greatest common depth 
between stations was used. Brunt-Váisálá (or buoyancy) frequency N was computed for 
each cruise at the stations nearest moorings P2 and P3. Usually the CTD cast was 
within a few kilometers of the moorings’ assigned position. Values for N were deter- 


mined by 


O) 


E l 
Neg SS] (18) 


© 


Values of p used in equation (18) were computed using the International Equation of 
State 1980 for seawater (Millero and Poisson, 1981) using S, T and P values from the 
CTD cast and were smoothed by center averaging over 10 m increments. Calculations 
of N? used Az = 10 m. The results from each CTD cast were plotted as a vertical profile 
of N? against depth. Most N? profiles revealed the main pycnocline at depths between 
50 m and 200 m. The highest value of N, 2.147 x 10? cps (a 7.76 minute period), was 
found near the main pycnocline. Minimum N values were found at the bottom of each 
cast with values near 2.344 x 10-* cps (a 71.01 minute period) and were assumed to be 
constant from the deepest calculated value to the sea floor. 

Profiles of the maximum trapping frequency w, = Nsina were constructed with the 


trapping frequency plotted in units of cpd against depth in meters. The possibility of 
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bottom trappmg at any particular CTD station 1s best indicated by the magnitude of the 
maximum trapping frequency calculated over the lowest 10 m of CTD data (Figure 10). 
The deepest calculated w,, typically 50 m off the bottom, was extrapolated to the bottom 
using the assumption that N? remained constant. Any wave with a frequency less than 
or equal to w, may be trapped. When @, is greater than 2.0 cpd, trapping may occur 
at both diurnal and semidiurnal frequencies. When @, 1s less than 1.0 cpd, no trapping 


of tidal currents should be expected. 


25 





AUN VA SACA REQUENA ROLE Nsina FAazQUENCY PROFILE 











| N**2 (10E-OE) CYCLES PER DAY 
0 5 10 16 20 25 30 35 40 46 60 o 6 10 16 20 
| O_o OOo o 
EE 
3 J Sl: 
3 | | 
dk? 
H + | 
28 se 5 
E Re 
| aa. n | 6 
E, Es ( 
EN SIE) 
58. EEN | 
E dH | 
E S. | | 
i T d | 
A et 
= | "i | 
E S 
= CN 
o i o | : 
A 8. ej: 
, = "x 
EZ 
| PROFILE 1D: CUmARSO STA 11 MOCRING P3 == 
se. 
É = | 
Val t 
Figure 10. Brunt-Váisálá Frequency and Trapping Frequency Profiles: 


N? and N sing are plotted versus depth using hydrographic data ob- 
tained near P3 in March 1990. 
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IV RESULTS 


pee CURRENT ELLIPSES 
Monthly averaged tidal current ellipses were constructed using the amplitudes and 
phases obtained from the complex demodulation of current meter records (Figures 11 
to 14). The plotted ellipses display a great deal of variability in their size, orientation 
and direction of rotation. The size and energy represented by the current ellipses was 
found to be greatest for both tidal frequencies at the deepest current meter of mooring 
p. 
l. Diurnal Current Ellipses 

Diurnal current ellipses for mooring P2 were constructed for seventeen consec- 
utive months (Figures 11 and 12). The energy represented by the ellipses was calculated 
using the sum of the squared axis lengths to compute kinetic energv in units of 
Em sec~ . At the 100 m level the average energy of the ellipses was 18.0 cm? sec?. The 
direction of rotation of 76% of these ellipses was counterclockwise. At 350 m, the av- 
erage energy dropped to 14.5 cm? sec~? with 70% of the ellipses rotating counterclock- 
wise. At 500 m, the average ellipse energy increased to 22.3 cm? sec”? with 65% rotating 
counterclockwise. The diurnal current ellipses were larger and more energetic at 500 m 
than at 350 m for each of the seventeen months analvzed and there were twelve months 
when more energy was present at 500 m than at 100 m. During July 1989 the diurnal 
current rotation at all three P2 current meters was clockwise and there were eight 
months when the rotation at three depths was uniformly counterclockwise. 

The current meter records for six months were used to construct diurnal current 
ellipses at mooring P3 (Figure 12). At the 100 m level the average energy represented 


by the ellipses was 25.1 cm? sec~? with 66% of the ellipses rotating counterclockwise. 
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Figure 11. Diurnal Current Ellipses from May 1989 to November 1989 at Mooring 
P2: Horizontal ellipses are shown with north up and the direction of 
rotation indicated by arrows. A velocity scale is indicated at the bot- 
tom. 


28 


O 


HIHON Q34VLN3IHO A X 


N 
S> 
NN NS 


07 
31V9S ALIDONIA Pd / 


06d38 065nYv oe1nr 08NNf OBÁVW O6udV OBUVIN 06833 O6NYTC 68930 


S3SdI 113 INI44NI IVN4NIG 


WOOOL 


WOOS 


INOSE 


NOO! 


WOON 


NOg€ 


WOOL 


Ed DNIBOOW 


td DNIHOOW 
Diurnal Current Ellipses from December 1989 to September 1990 at 


Moorings P2 and P3: Details are as in Figure 11. 


Figure 12. 
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Figure 13. Semidiurnal Current Ellipses from May 1989 to November 1989 at 
Moorings P2 and P3: Details are as in Figure 11. 
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Semidiurnal Current Ellipses from December 1989 to September 1990 


at Mooring P2 and P3: Details are as in Figure 11. 


Figure 14. 
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At 350 m, the average energy dropped to 12.4 cm? sec? . At 500 m, four months were 
analyzed and resulted in an average energy of 13.5 cm?sec? with 75% rotating 
counterclockwise. The smallest ellipses and lowest energy (9.5 cm? sec?) were found at 
1000 m with 66% of the ellipses rotating counterclockwise. The ellipses were larger and 
contained more energy at 500 m than did those at either 350 m or 1000 m. There were 
four months when the current rotation at each P3 current meter was counterclockwise. 

The bathymetric contours at mooring P2 are oriented at approximately 005 *T 
and the average orientation of the 500 m diurnal tidal current ellipses at P2 was 049* 
T. The semi-major axis of the diurnal ellipses crossed the isobaths a mean relative angle 
of 44°, At mooring P3 the bathymetry is oriented at 020° T. The average orientation 
of diurnal tidal ellipses at 1000 m was 074° T and the semi-major axis of these ellipses 
crossed the isobaths at a mean relative angle of 54°. 

2. Semidiurnal Current Ellipses. 

The semidiurnal current ellipses at P2 and P3 were larger and more energetic 
than those at diurnal frequencies. The orientation of the semidiurnal ellipses was less 
variable than the diurnal ellipses with almost all of them oriented in a north-northeast 
direction. 

At mooring P2, the seventeen month average energy represented by the ellipses 
at 100 m was 59.1 cm? sec~*. The rotation of 58% of these ellipses was counterclockwise. 
At 350 m, the average energy decreased to 48.2 cm? sec”? and 52% of the currents ro- 
tated counterclockwise. At 500 m the average energy increased to 85.4 cm? sec”? with 
47% of the ellipses displaying counterclockwise rotation. The semidiurnal currents dis- 
plaved the same trends exhibited by the diurnal current ellipses where the largest ellipses 
and greatest energies were found at the deepest level and the smallest ellipses and least 


energy was at 350 m. The rotation of semidiurnal current ellipses at P2 was clockwise 
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at all three current meters during four months and counterclockwise during four other 
months. The remaining nine months contained ellipses with opposing rotation. 

The semidiurnal current ellipses at mooring P3 were notably smaller than those 
at P2. The average energy represented by the ellipses at 100 m was 28.4 cm? sec”? with 
66% rotating clockwise. At 350 m the average energy decreased to 22.7 cm? sec? with 
66% of the currents rotating counterclockwise. A four month average at 500 m resulted 
in an average semidiurnal energy of 34.5 cm? sec~? with all four ellipses rotating clock- 
wise. At 1000 m the average energy was only 8.3 cm? sec. At this level, 66% of the 
ellipses rotated clockwise. Like the diurnal currents at P3, the semidiurnal currents were 
greatest at 500 m. The smallest ellipses and least energy was at the 1000 m level. At 
mooring P3 two of the six months analyzed had semidiurnal currents at four depths ro- 
tating in a clockwise direction and there was one month when the rotation was 
counterclockwise at each depth. 

The average orientation of 500 m semidiurnal current ellipses at P2 was 030° T 
and the isobaths were crossed at an angle of 25°. The average orientation of 1000 m 
semidiurnal current ellipses at P3 was 043° T and the semi-major axis crossed the 


isobaths at 23°. 


B. ENERGY PROFILES 

Vertical profiles of energy at both the diurnal and semidiurnal frequencies (Figures 
15 and 16) were constructed for each month using the output from the spectral analysis 
program. These plots contain information similar to that of the tidal current ellipses 
expressed in a different format. At P2 the highest energy levels were at the deepest 
current meter for both frequencies. The energy levels at P3 show low energy levels at 


the deeper current meters. 


33 


NOV 89 


OCT 89 














O) 

oo 

a 

Lu 

LO I 
« 
z 
ecc 

o D 

co o 

9 E 
N 

< 

eo 

cO 

E 

2 

5 

O) 

: z 

Z a 

2 2 

73 a 





MAY 89 





Figure 15. Spectral Energy Profiles from May 1989 to November 1989 at Mooring 
P2: All plots are scaled identically. 
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Figure 16. Spectral Energy Profiles from December 1989 to September 1990 at 
Moorings P2 and P3: All plots are scaled identically. 


35 


1. Diurnal Energy Profiles 

At mooring P2 the average monthly diurnal energy levels were 13.9 cn? sec? at 
100 m, 13.4 cm? sec”? at 350 m and 20.2 cm? sec”? at 500 m. The energv at 500 m ex- 
ceeded that at 350 m for each of the seventeen months analvzed except September 1989. 
The energy at 500 m was also higher than that at 100 m except for the months of June 
1989, February 1990 and April 1990. The curves plotted in Figures 15 and 16 exhibit a 
‘C’ shape with minimum energy levels tvpically found at a depth of 350 m. 

Diurnal energv at P3 was greatest near the surface. Average monthly values 
were 17.3 cm? sec”? at 100 m, 7.2 cm? sec? at 350 m 9.7 em sec? at 500 m Ed 
cm? sec”? at 1000 m. The curves plotted below 100 m are nearly uniform and show vir- 
tually no amplification of energy toward the bottom. 

2. Semidiurnal Energy Profiles 

The semidiurnal current energy levels at P2 were likewise greatest at the deepest 
current meters. The same 'C' shaped curve observed with diurnal currents, with a min- 
imum energy at 350 m, Was exhibited by the semidiurnal currents. Average monthly 
values of 63.1, 47.6, and 77.4 cnv sec? were found at 100 m, 350 m and 500 m respec- 
tively. The energv at 500 m was greater than that at 350 m for fourteen of the seventeen 
months analyzed. The energv at 500 m exceeded that at 100 m for nine of the seventeen 
months. 

At P3, the average semidiurnal energy levels were fairly consistent until drop- 
ping off sharply at 1000 m. Values of 20.1, 21.3, 30.3 and 7.3 cnv sec? were found at 
100 m, 350 m, 500 m and 1000 m. No bottom amplification of semidiurnal energy was 


evident at P3. 


C. TEMPORAL VARIABILITY 
The same data used to plot the vertical profiles were used to create time series 


showing the change of energy at each depth for the two moorings. Diurnal bottom 
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trapping was evident at P2, with energy at 500 m greater than the energv at 350 m for 
the entire seventeen month period analvzed (Figure 17). Bottom trapping Was always 
present and was most pronounced during July - August 1989, October 1989 - January 
1990 and Mav - September 1990 when the energy at 500 m exceeded that of any other 
depth. The magnitude of the diurnal energy at 500 m was statistically greater than 350 
m or 100 m at 95 percent confidence for the period October 1989 through January 1990. 
At mooring P3 no diurnal bottom trapping was indicated by the energy at 350 m, 500 
m or 1000 m (Figure 18). Diurnal energy did not increase toward the bottom and the 
energy near the surface, at 100 m, was alwavs greater than at deeper current meters. 
Bottom trapping of semidiurnal energy occurred at P2 for fourteen of seventeen 
months when the energv at 500 m exceeded the energy at 350 m (Figure 19). In addition, 
the energy at 500 m was greater than 100m for nine of the seventeen months analyzed. 
Bottom trapping at 500 m was statistically significant using a 95 percent confidence in- 
terval for September 1989, December 1989, April - May 1990 and September 1990. At 
P3 no semidiurnal trapping was evident and the energy at 1000 m was less than at any 


other depth for five of the six months analyzed (Figure 20). 


D. APPLICATION OF TRAPPING THEOR Y 

The size of the current ellipses and energy distribution at each depth suggest that 
an amplification of energy occurs near the botton of mooring P2. A similar pattern Was 
not obvious at P3. To compare observed motions to theoretical trapped waves an ex- 
ponential curve, representing the decay of trapped energv, was compared to the observed 
profile of kinetic energy. The theoretical curves were constructed using equation (12) 
with x given by equation (16). 

]. Trapping Scales and Decay Curves 

At mooring P2 the maximum trapping frequency is exactly equal to the diurnal 


frequency (w, = 1 cpd) when, by equation (7), N = 1.930 x 10º cps . Using H = 800 m 
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Figure 17. Time Series of Diurnal Energy at Mooring P2: The temporal vari- 
ability of kinetic energy is shown with curves representing the energy 
levels at current meter depths. 
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Figure 19. 
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and f=1.375x 10-5 cps in equation (16) results in a trapping coefficient of 


— 3.3116 x 10? mr'. The horizontal length scale of this wave is L = + = 26.632 km. 


Similar calculations were conducted for the first five baroclinic modes at both P2 and 
P3 for diurnal and semidiurnal trapping using N values determined by equation (7). The 


coefficients are tabulated in Tables 3 and 4. 


Table 3. TRAPPING COEFFICIENTS AT MOORING P2 












Table 4. TRAPPING COEFFICIENTS AT MOORING P3 
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For the baroclinic modes, the trapping coefficients in Tables 3 and 4 were used 
to produce a vertical trapping scale. The trapping scale (or e-folding scale) indicates the 
height off the bottom where the energy of a trapped wave decayed by a factor of e! . 
Most of the trapped energy for each wave was contained in a bottom layer with a 
thickness (determined by the trapping scale) indicated in Table 5. Note that the two 


deepest current meters at P2 (located 450 m and 300 m from the bottom) fell within this 
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laver but the lowest instrument at P3 (located 800 m from the bottom) did not. For 
higher modes trapping 1s stronger, trapping scales are smaller and energy decays more 


rapidly from the bottom. 


Table 5. TRAPPING SCALES AT NOORINGS P2 AND P3 






Using the appropriate trapping coefficient, exponential decay curves for each 
mooring were constructed which represent the first two baroclinic modes at diurnal and 
semidiurnal frequencies at moorings P2 and P3 (Figures 21 and 22). Similar curves were 
fit to observed spectral energy profiles to determine if the vertical energy distribution off 
Point Sur indicated the presence of bottom trapped waves. The curves were constructed 
using N values determined from CTD data whenever the diurnal or semidiurnal fre- 
quencies were less than the maximum trapping frequency. 

2. Analysis of Bottom Trapping 

CTD data from each of the hydrographic cruises in Table 2 was analyzed to see 
if bottom trapping was possible at P2 or P3 at diurnal or semidiurnal frequencies using 
the criteria that w, = N sina must be equal to or greater than the tidal frequencies. 
Each CTD cast taken at a mooring site was assumed to be representative of the aver- 
aging interval used to compute the spectral estimates (one month). 

When bottom trapping was possible, theoretical decay curves were constructed 
using observed values of H, N, and appropriate trapping coefficients (x). The curves 
were compared to the observed energy levels obtained from spectral analysis by match- 
ing the curves at the level of the deepest current meter (500 m at P2) and comparing the 


theoretical and observed energy at other (100 m and 350 m at P2) levels. When the 
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Figure 21. Theoretical Trapping Profiles at Mooring P2: The exponential decay 
of trapped energy (with arbitrary units) is illustrated for the first and 
second modes at tidal frequencies. 


modeled energy fell within the 95 percent confidence intervals of the observations, the 
change of energy with depth was considered to be consistent with bottom trapped waves. 
The buoyancy frequencies N, maximum trapping frequencies w, and trapping 


coefficients x, for each month when data were available are tabulated in Table 6. The 
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Figure 22. Theoretical Trapping Profiles at Mooring P3: The exponential decay 


of trapped energy (with arbitrary units) is illustrated for the first and 
second modes at tidal frequencies. 


maximum trapping frequency w, at P2 exceeded the diurnal frequency (w = 1.0 cpd) 
each of the eight months. The maximum trapping frequency exceeded the semidiurnal 


frequency (w = 2.0 cpd) except for the month of January 1990 (w, = 1.22 cpd). 
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Table 6. OBSERVED BUOYANCY FREQUENCIES, TRAPPING FREQUEN- 
CIES AND MODE I! TRAPPING COEFFICIENTS AT MOORING P2 












The diurnal energy distribution at P2 (Figures 23 and 24) agreed with the mod- 
eled decay for seven of the eight months. During January 1990, the observed energy at 
350 m was less than indicated by theory. The theoretical curves, extrapolated to a depth 
of 800 m, indicated an average ee level of 40.5 cm? sec”? trapped at the sea 
floor. 

The maximum trapping frequency w, at P2 exceeded the semidiurnal frequency 
(c, — 2.0 cpd) during seven of the eight months. The observed energy at 100 m during 
May 1989 was higher than was modeled by matching the theoretical curve to the 500 
m observation. The modeled curve fell within the 95 percent confidence intervals for the 
six remaining months. Extrapolation of the theoretical curve to 800 meters produced 
an average kinetic energy level of 310.0 cm? sec”? at the bottom over the time period 
studied. 

The buoyancy frequency calculated over the deepest 10 m of CTD data ranged 
from 20.3 cpd to 66.2 cpd and the magnitude of the maximum trapping frequency, which 


was proportional to N , ranged from 1.22 cpd to 3.97 cpd. The fluctuations of N had a 


46 


much smaller effect on the trapping coefficients. The diurnal trapping coefficient at P2 
varied only 9.0 x 10-* m*'.. The semidiurnal trapping coefficient had a range of only 
2.1 x 107 m-'. The consistent trapping coefficients produced nearly identical trapping 
scales at P2. The diurnal trapping scale (defined by the e-folding scale) was 302 m and 
the semidiurnal trapping scale was 151 m for each of the eight months analyzed. 

The diurnal energy distribution at P3 did not have a profile indicative of bottom 
trapped waves. Since the maximum trapping scale for a diurnal wave at P3 was ap- 
proximately 151 m and the deepest current meter at P3 was 800 m above the sea floor, 
any energy from a trapped wave would have decayed 99.5 percent before reaching this 
current meter and could not be observed. The maximum trapping frequency at this 


mooring never exceeded w, = 2.0 cpd and no semidiurnal trapping occurred. 
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Figure 23. 





Trapping Frequencies and Energy Profiles at P2 During 1989: Maxi- 
mum trapping frequencies are plotted against depth with the dashed 
line indicating 2.0 cpd in the top row. The middle row shows observed 
diurnal spectral energy profiles (solid lines) with 95% confidence inter- 
vals (horizontal bars). Theoretical energy profiles (dashed lines), based 
on a decav scale calculated using N from the bottom of CTD casts (or 
at $00 m for Julv 1989) , are matched to the observations at 500 m. 
Sinular profiles are shown in the third row for semidiurnal energy. 
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Figure 24. Trapping Frequencies and Energy Profiles at P2 During 1990: Data is 
plotted as described in Figure 23. No theoretical semidiurnal energy 
profile is shown for January 1990 when N sina near the bottom was 
less than 2.0 cpd. 
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V. DISCUSSION 


A. BOTIOM TRAPPED WAVES 

The trapping or amplification of tidal band kinetic energy near the ocean bottom 
was evident from the data obtained at mooring P2. Both the current ellipses and the 
energy distribution displaved a significant increase from a depth of 350 m to 500 m 
during the seventeen months of data analyzed. A simular increase was not evident from 
the data obtained at mooring P3. 

The two methods of analysis used to examine the currents produced results which 
were in general agreement. The energy levels derived through complex demodulation 
and current ellipses Were larger in magnitude than those obtained through spectral 
analysis. The total KE represented by the ellipses was 10.7 % greater than the energy 
obtained from spectral analysis. The lower energy spectral analysis values may be due 
to Doppler shifting of small wavelength internal tides to other parts of the spectrum. 
This shifting could spread energy around the tidal band so that energy measured pre- 
cisely at the frequency of a tide, using the high resolution FFT of VCSPC3, improperly 
excluded tidal energy located in sidebands (W’unsh, 1975). The difference between the 
energy from ellipses and spectral analysis was greatest at P3 during December 1989. The 
spectral density plots for the current meter records at 100 m and 350 m show that the 
diurnal peak was spread out over a wider frequency band (perhaps by Doppler shifting) 
and had a smaller amphtude than for other months of the same record. 

Both ellipses and spectral energy profiles indicated more energy was present at 100 
m than at 350 m at both moorings for diurnal and semidiurnal currents. The current 
meters at 100 m were above the main pycnochne and the higher energy levels in the 


surface laver could be the result of wind forcing at or near tidal frequencies. Strong 
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coastal winds with diurnal variations have often been observed along the California 
coast (Beardsley et al., 1987). The pycnocline could prevent coupling of surface energy 
to deeper layers. Diurnal tidal currents may also be influenced by inertial currents with 
a frequency is 1.1 cpd off Point Sur. The narrow frequency separation between diurnal 
and inertial currents could result in smearing or leakage of energy into the diurnal band. 

The observed distribution of tidal energy had a vertical decay rate indicative of 
bottom trapping at mooring P2. The fit of a theoretical curve to the data fell within the 
95 percent confidence intervals with few exceptions. In January 1990 the observed 
diurnal energv at 350 m was less than predicted by bottom trapped theory using the 
buoyancy frequency calculated from a single CTD cast. If the buoyancy frequency had 
been greater the trapping coefficient, x, would have been larger. A larger k would result 
in a faster decay rate and the theoretical energy at 350 m would have been closer to 
observations. The semidiurnal energy at 100 m during Mav 1989 was higher than ex- 
pected using the modeled decay curve and these energetic semidiurnal currents may re- 
sult from surface (wind) forcing. Another factor in the vertical semidiurnal energy 
distribution is free internal waves which exist at superinertial frequencies. Free internal 
waves could be generated by tidal forcing across the continental shelf break and could 
then Eus into deeper stratified water (Baines, 1982). Horizontal velocities over the 
shelf have been observed where the highest velocities were near the surface and bottom. 
The velocity structure was attributed to the first dynamical mode (Rosenfeld, 1990) and 
the same structure may dominate internal tides over the continental slope. The 
semidiurnal currents at P2, located close to the shelf break, were twice as energetic as the 
semidiurnal currents found further offshore at P3 and did have more energy at the top 
and bottom current meters (100 m and 500 m) than at the middle (350 m). Bottom 


trapping did not occur with anv discernable seasonal pattern. The absence of distinct 
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temporal variation agrees with other observations of internal tides along the California 
coast (.Voble et al., 1987). 

The theoretical curves plotted for P2 in Figures 23 and 24 were extrapolated to 800 
m. The average diurnal energy represented by these estimates was 40.5 cm? sec-2. Using 
the average diurnal current ellipse eccentricity (e = 3.42), diurnal currents at 800 m 
would have had a semi-major axis of 6.] cm sec! and a semi-minor axis of 
1.8 cm sec-!. The theoretical energy for semidiurnal currents extrapolated to 800 m was 
310.0 cm? sec?. An average semidiurnal eccentricity of e = 3.26 produced an ellipse re- 
presenting bottom currents with semi-major and semi-minor axis lengths of 16.8 
cm sec”! and 3.2 cm sec”! respectively. The magnitude of these currents is reasonable, 
however, bottom trapped theory assumes inviscid boundary conditions and bottom 
friction may reduce the velocity of predicted bottom trapped currents. 

The average wavelength of bottom trapped diurnal waves, computed using the av- 
erage trapping coefficient in Table 5, was 66.5 km. The average wavelength of bottom 
trapped semidiurnal waves was 33.2 km. This compares favorably to first baroclinic 
mode internal waves of semidiurnal frequency observed during CODE which had wave- 
lengths between 20 and 30 km (Rosenfeld, 1990). 

The rotation of tidal current ellipses is not described by bottom trapped theory. The 
observed rotation of the current ellipses was variable with 73 percent of the diurnal el- 
lipses rotating counterclockwise and 52 percent of the semidiurnal ellipses rotating 
clockwise. Barotropic Kelvin waves, free internal waves and shelf waves all influence the 
rotation of current ellipses at tidal frequencies. The observed direction of rotation 1s 
determined by the cumulative effect of these unresolved wave regimes. 

Barotropic semidiurnal tides are dominated by the M, constituent. The horizontal 
velocities generated by Kelvin waves at this frequency are normally aligned with topog- 


raphy (Noble et al., 1987). The observed tidal current ellipses all had cross-slope com- 
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ponents and vertical velocity structures which may reflect the influence of a baroclinic 
wave structure. At the deepest current meters, the diurnal current ellipses were oriented 
more across the bottom slope (P2 = 44°, P3 = 54°) than were the semidiurnal ellipses ( 
P2 = 25°, P3=23° ). With this geometry the particle motion of diurnal currents had a 
larger upslope component and more bottom trapping than semidiurnal currents of the 
same magnitude. Even though observed diurnal currents were significantly less energetic 
than semidiurnal currents, diurnal bottom trapping was always evident at mooring P2. 

Both moorings had a similar bottom slope « and the currents analyzed were limited 
to diurnal and semidiurnal frequencies. The only variable in determining the limiting 
trapping frequency ( w,= N sina ) was the buoyancy frequency. Higher values of N 
resulted in higher trapping frequencies. For each of the eight months when CTD casts 
were made at P2, diurnal trapping was possible as indicated by the calculated values of 
©, and semidiurnal trapping Was indicated during seven of these months. The current 
ellipses and energy profiles exhibited bottom trapping in each case. No bottom trapping 
was evident at P3 although there were two months with o, values high enough to permit 
the trapping of diurnal energy. The absence of any notable bottom trapping at P3 may 
be due E the trapping scale and the depth of the instruments at that mooring. All three 
current meters at P3 were above the e-folding depths of bottom trapped diurnal and 
semidiurnal energy. If any trapped energy was present, the magnitude of the signal was 
indiscernible at the P3 current meters. 

The difference in maximum trapping frequencies at each mooring is caused by dif- 
ferent hydrographic conditions at P2 and P3. The cruise of November 1989 can be used 
to illustrate typical conditions. Density gradients are greater near the ocean bottom at 
P2 than at P3 (Figure 25) and the gradients produce higher N values at P2 than those 
calculated at P3. The stronger gradients at P2 permit trapping of waves with frequencies 


higher than the maximum frequency of trapped waves at mooring P3. 
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Figure 25. Vertical Density Section Off Point Sur, California During November, 
1989: The positions of moorings P2 and P3 are indicated. 


D. LIMITATIONS 

Several assumptions were made in comparing the observations at P2 with theoretical 
predictions for bottom trapped waves. The first was that the tidal currents were always 
oriented across the sea floor slope so that the angle ¢ = 90°. This angle permits maxi- 
mum trapping and simplifies the equation o = N sina sino to w,=N sina. At angles 
less than 90°, the limiting trapping frequency c, will be lower and the bottom trapping 


of tidal energy may not occur. lf was assumed to be the angle at which the semi-major 
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axis of the tidal current ellipses crossed the bottom contours, then the maximum trap- 
ping frequency c at P2 ($ 2 44* at 500 m) would still exceed 1.0 cpd for seven of the 
eight months when N values were available. For the P2 semidiurnal ellipses ( 
p = 25” at 500 m), o Would not have exceeded 2.0 cpd for any month. If no semidiurnal 
trapping was possible, then the vertical distribution of energy at semidiurnal frequencies 
must be attributed to other factors including free internal waves. 

A second assumption was that a single CTD cast provided information which could 
be used to interpret monthly averaged current meter data. The CID casts typically 
ended 50 m above the sea floor and salinity and temperature values over the deepest 10 
m interval were used to calculate N for each month. These N values were assumed to 
be constant from the bottom of the CTD cast to the sea floor. Trapping frequencies 
calculated with N were similarly extrapolated. The use of extrapolated values was sup- 
ported by the CTD cast of July 1989 which was slightly offshore of mooring P2 and ex- 
tended to 903 m (the 10 m interval centered at 800 m was used to calculate N at P2). 
The values of N were fairly constant during July 1989 (42.2 cpd - 40.9cpd) at the depths 
where extrapolation was used (650 m to 800 m) for the seven other CTD casts. As dis- 
cussed in the previous section, variations in N could affect the fit of the theoretical decay 
curves to observations. The c, and k values using a constant N were compared to 
©, and k computed using N obtained by extending the slope of the N? curves to 800 m. 
The differences created slight variations of œ, which did not conflict with any of the re- 
sults presented in Chapter IV. The differences in the x values were very small and con- 
sidered insignificant. 

The instrumentation at each of the moorings was not ideally positioned for the study 
of bottom trapped waves. The current meters were too sparse to provide a good profile 
of the currents throughout the water column. The current meters at the 100 m level were 


above the pycnocline and may have contained energy from surface forcing which was 
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not considered in constructing the theoretical decay curves. The lowest current meters 
at P2 and P3 were not deep enough to detect the strongest currents expected near the 
bottom. At P3 the small trapping scales made it unlikelv that any signal from a bottom 
trapped wave could be detected. The current meter at 1000 m was well above the e- 
folding depth of diurnal and semidiurnal waves. Additionally, the location of P3 in an 
area with relatively weak density gradients restricted most bottom trapping to subtidal 


frequencies. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


The analysis of current meter records from moorings located off Point Sur, 
California showed an increase of energy at tidal frequencies with depth. This can po- 
tentiallv be explained using the theory of Rhines (1970) OO trapped waves over 
the continental slope. The increased energv was observed whenever bottom slope, water 
stratification and frequency met prescribed conditions. 

Seventeen months of current meter data was examined from mooring P2. The data 
were examined using complex demodulation and spectral analvsis. At mooring P2, 
bottom trapped energy at the diurnal frequency was observed for each month analyzed. 
Bottom trapping occurred at semidiurnal frequencies for fourteen of the seventeen 
months but the energy contribution due to free internal waves, which influenced the 
semidiurnal results, Was not determined. Bottom trapping at mooring P2 occurred dur- 
ing each month of the year and no seasonal patterns could be distinguished. The maxi- 
mum trapping frequency at P2, computed using N sina, was an accurate indicator of 
bottom trapping at tidal frequencies. 

No bottom trapping was observed at mooring P3 using data collected over six 
months. Weak density gradients near the sea floor often kept the maximum trapped 
frequency below the tidal frequencies. The current meters at this mooring were posi- 
tioned too far above the sea floor for detection of bottom trapped waves at diurnal or 
semidiurnal frequencies. 

The bottom trapped wave theory described by Rhines (1970) was also used to cal- 
culate vertical decay scales and energy distribution profiles. The theoretical decay was 


modeled using hydrographic data collected near the mooring locations. In most cases 


Dy 


the modeled decav agreed with the observed decay of tidal energy away from the sea 
floor within the error of the calculation. 

The tidal currents off Point Sur are the most energetic currents found over the con- 
tinental slope. The semidiurnal currents dominate but diurnal currents can be signif- 
icant, particularly when amplification due to bottom trapping occurs. Further study of 
bottom trapping should be conducted using the continuing current meter arrays at or 
near the location of P2. Additional deep instrumentation, particularly just above the sea 
floor, is needed to adeguately study this phenomena. Because trapping Was observed 
throughout the year, a short deplovment of a few months length may collect sufficient 
data to verifv trapping at tidal frequencies. Tidal currents are at the upper limit of the 
range of trapped waves and longer deplovments could be used to examine coastally 


trapped waves occurring at much lower frequencies under the same conditions. 
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